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ABSTRACT 



The four year X-ray all-sky survey (eRASS) of eROSITA telescope aboard the Spektrum-Roentgen-Gamma satellite will detect ~ 3 
million AGN with a median redshift of z « 1 and typical luminosity of L0.5-2.0kev ~ 10 44 erg s -1 . We show that this unprecedented 
AGN sample, complemented with redshift information, will supply us with outstanding opportunities for large-scale structure re- 
search. For the first time, detailed redshift and luminosity resolved studies of the bias factor for X-ray selected AGN will become 
possible. The eRASS AGN sample will not only improve the redshift and luminosity resolution of these studies but will also expand 
their luminosity range beyond L0.5_2.0kcV ~ 10 44 erg s -1 , thus making possible direct comparison of clustering properties of lumi- 
nous X-ray AGN and optical quasars. These studies will dramatically improve our understanding of AGN environment, triggering 
mechanisms, growth of super-massive black holes and their co-evolution with dark matter halos. 

The eROSITA AGN sample will become a powerful cosmological probe. It will make possible detection of baryonic acoustic oscilla- 
tions (BAOs) for the first time with X-ray selected AGN. With the data from entire extragalactic sky, BAO will be detected at a S 10<x 
confidence level in the full redshift range and with ~ 8<x confidence in the 0.8 < z < 2.0 range, currently uncovered by any existing 
BAO surveys. In order to exploit the full potential of the eRASS AGN sample, photometric and spectroscopic surveys of large area 
and a sufficient depth will be needed. 

Key words. Surveys - X-rays: general - Quasars: general - Galaxies: active - Cosmology: Large-Scale Structure of Universe 



1. Introduction 

Large-scale structure (LSS) studies are established as an impor- 
tant tool for studies in two major areas of astrophysics: cos- 
mology and galaxy evolution. A key of their success is the in- 
creasing number of surveys at different wavelengths with in- 
creasing depths and sky coverages. In X-rays, many deep, ex- 
tragalactic surveys have been performed in the last decade 
(Bran dt & Hasingerl I2005I) . However, in comparison to other 
wavelengths, X-ray surveys with a large sky coverage and with 
sufficient depth are still rare. Th e previous X-ray all- sky sur- 
vey w as performed by ROSAlQ (Truemperl ll993t IVoges etaTI 
1999) more than two decades ago. Its successor with a ~ 30 
times better sensitivity will be the four year long all-sky survey 
(eRASS) of the eROSITAEI telescope (PredehTetaDIIOlO), to be 
launched aboard the Russian Spektrum-Roentgen-Gamma satel- 
lite in 2014. 

The major science goals of the eROSITA mission are to 
study cosmology with clusters of galaxies and active galactic nu- 
clei (AGN) and to constrain the nature of dark matter (DM) and 
dark energy. For a comprehensive descript ion of the eROSIT A 
mission see the science book of eROSITA dMerloni et al.l2012r) . 

In this work we explore the potential of studying LSS with 
the AGN sample to be detected in eRASS. We focus on two im- 
portant aspects of LSS studies: the clustering strength (repre- 
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sented by the linear bias factor, Sect. |3) and the baryonic acous- 
tic oscillations (BAOs, Sect. [4]!. For measuring the former quan- 
tity redshift accuracy of photometric surveys is sufficient, there- 
fore bias studies can be successfully conducted during and soon 
after the time eRASS is finished. The BAO measurements, on 
the other hand, will be much more difficult to accomplish as the 
spectroscopic redshift accuracy over large sky areas will be re- 
quired. Note that sufficient redshift accuracy can be also pro- 
vided by high quality narrow-band multi-filter photometric sur- 
veys. 

In our previous work dKolodzig etal.ll2012h . we have studied 
the statistical properties of the AGN sample of eRASS and will 
adopt these results here. In the current work, we focus on the 
AGN detected in the soft energy band (0.5 - 2.0 keV) and on 
the extragalactic sky (\b\ > 10°, ~ 34 100 deg 2 ). In the following 
calculations we will assume t he four year average se nsitivity of 
1.1 x 10~ 14 erg s 1 cirT 2 from lKolodzig et all d2012l Table 1). 

Optical followup surveys will be needed to provide identi- 
fication and redshift information to a desired accuracy for all 
eRASS AGN. Current optical surveys are not sufficient in size 
and/or depth. A sensitivity of / ~ 22.5 mag (R ~ 23.0 mag) is re- 
quire d to detect at least 95 % of the eRASS AGN ( Kolodzi g~eTaT] 
2012). Many photometric and spectroscopic surveys with dif- 
ferent parameters have been proposed or are being already in 
construction (e.g. lMerloni et al.l (120121) ). For the purpose of our 
investigation we will assume that redshifts are available for all 
eRASS AGN. We will explore the effects of red shift-errors in 
the forthcoming paper (lHutsietal.112013. in prep.). 
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Fig. 1. The angular power spectrum of the full eRASS AGN 
sample (soft band, four years) for the extragalactic sky (/ s k y ~ 
0.83) and < z < 5. The gray shaded area and the blue his- 
togram show the lcr uncertainty region (Eq.|4) without and with 
/"-binning, respectively. The horizontal dotted line shows the 
level of shot noise, which was already subtracted from the an- 
gular power spectrum. For multipoles above the vertical dashed 
line (representing Z max * 500) our assumption of a linear clus- 
tering starts to break down. Therefore, we do not consider these 
multipoles in our further calculations. 
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Fig. 2. Same as Fig. Q]but additionally with the angular power 
spectra for various narrow and broad redshift ranges (see 
Sect.l2T2l. 



We assume for this work a flat ACDM cosmology with the 
following parameters: Hq = 70kms~' Mpc~' (h = 0.70), Q m = 
0.30 (Qa = 0.70), Qb = 0.05, <x 8 = 0.8. Luminosities are given 
for the soft energy band (0.5 - 2.0 keV) and we use the decimal 
logarithm throughout the paper. 



2. Angular power spectrum 

The common tool to study LSS is to measure and analyze the 
clustering of objects (such as AGN) with the 2-point correla- 
tion f unction (2pCF) or the power spectrum (PS) (e.g. [Peebles! 
1980). The two me thods, 2pCF and PS, h ave their benefits and 
disadvantages (e.g. Wa ll & Jenkinsll2"oT2h but should contain in 
the end the same information about the LSS. In this work we 
will use the angular power spectrum Q to characterize cluster- 
ing properties of objects. In order to predict the power spectra 
which will be measured with eRASS AGN, we will rely on the 



model for AGN clustering of iHiitsi et al.1 (1201 2l) . in particular, 
we will use their model II. The details of our calculations are 
summarized below. 

We calculate the angular power spectrum as follows: 



Q = ^ J P{k) [W ( (k)f k 2 dk 
where the projection kernel is 



W e (k) 



f(z) 8(z)b(M en ;z)j ( (kr(z))dz 



(D 



(2) 



Hereby, P(k) is the z — 3D li near power spect r um, fo r which 
we will use fitting formulae of Eisens tein & Hul (1 19981) . f(z) is 
the normalized radial selection function, g(z) is the linear growth 
function (e.g. Dodelson 2003), b(M eS , z) is the AGN linear clus- 
tering bias factor and jc are the spherical Bessel functions of 
order (, wh ere r(z) is the co-moving distance to redshift z (e.g. 
iHoggll 19991) . 

The radial selection function is defined as the (normalized) 
differential redshift distribution of AGN, which we calculate 
with the AGN X-ra y luminosity function (XLF), 4>(\og L, z), of 
lHasinger et al.l (20050. It is the only quantity that contains the 
information about eRASS, since it depends on the survey sensi- 
tivity [S ) as following: 



<W dV(z) 
— (S,z) = — — 
dz dz 



logL„,, 

/ 



0(logL,z)dlogL (3) 



Here, [Mpc 3 sr 1 ] is the co-moving volume element and 

L m i„ (S,z) = 4nS d 2 (z), where d^(z) is the luminosity distance 
re.g. lHog3[T999h . 

The AGN linear clustering b ias factor, b(M e tf, z), is com- 
puted with the analytical model of lSheth et al.l (12001) by assum- 
ing an effective mass M e jf of the DM halo (D MH) where the 
AGN reside. Based on re cent observations (e.g. Allevato et all 
l2^TTtlKrumpe et al.ll20ll iMountrichas et al]|2013lh we assume 
an effective mass of M eff = 2 x 10 13 h~ l Mq. 

In our work, we only focus on the linear clustering regime. 
Therefore we restrict our calculations to the spatial co-moving 
scales larger than k m . dx « 0.2 h Mpc -1 , corresponding to wave- 
lengths larger than ~ 30 hr l Mpc. The associated multipole num- 
ber is /"max© = &max f(z) an d depends on the median redshift z 
of the considered redshift bin. At the median redshift of eRASS 
AGN sample this is £ max (z ~ 1) ~ 500. Thus, for our calculations 
we do not consider C( at multipoles higher than ( mdx . 

For simpli city, we do n ot take linear redshift space distor- 
tions (RSD) (lKaiseilll987h into account. Since the signal-to- 
noise ratio (S/N) in our angular power spectra is rather poor at 
small multipoles (see Figs.Q]and[2]i, where linear RSD become 
most significant, we do not expect that our results would change 
significantly if we would consider them in our calculations. 



2.1. Uncertainties 

The variance of the Q can be well approximated with: 



(21 + 1) / sky 



(4) 



See lKolodzig et alJ (2012) for details. 
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assuming Gaussian statistics of the matter fluctuations (£ ^ 
Anax)- Hereby, is the sky fraction, which takes into account 
the effective loss of modes due to partial sky coverage, and N 
is the AGN surface number density [sr~'], which is computed 
with the AGN XLF and the survey sensitivity of eRASS (See 
iKolodzig et alj|20l"2l) . The first term (Q) in the brackets repre- 
sents the cosmic variance and becomes important at large scales 
(small (). The second term, the shot noise (Af -1 ), takes into ac- 
count that we are using a discrete tracer (AGN) and becomes 
dominant at small scales (large (), where N~ l » Q (see e.g. 
Fig- ID- In order to minimize the uncertainty in Q both high sky 
coverage and large number density of objects are needed. 



2.2. Results 

In Fig. Q] we show the expected angular power spectrum of the 
full eRASS AGN sample after four years for the entire extra- 
galactic sky. By introducing redshift information (available from 
other surveys), the angular power spectrum becomes a relevant 
tool for LSS studies. In Fig. [2] we can see how its amplitude 
increases with a decreasing size of the redshift bin and also os- 
cillations (see Sect. |4| in the angular power spectrum become 
more prominent. We can see from the two angular power spectra 
of 0.80 < z < 0.85 and 2.00 < z < 2.05 (with same redshift bin 
size) in Fig. [2] that the turnover of the spectrum and the posi- 
tions of the oscillations depend on the redshift. The amplitudes 
are also different because the clustering strength increases with 
redshift (see Sect. |3). Since the redshift distribution of eRASS 
peaks around z ~ 0.8, the number density around z ~ 0.8 is 
much higher than at z ~ 2.0 and therefore the uncertainty of the 
angular power spectrum is smaller for 0.80 < z < 0.85 than for 
2.00 < z < 2.05. 



3. Linear bias factor 

The linear bias factor b is an important parameter for cluster- 
ing analysis of AGN. It connects the underlying DM distribu- 
tion with the AGN population. Observationally, it has been so 
far a very challenging task to measu re this connection wi t h high 
accuracy, due to low statistics (e.g. iKrumpe et"aD 120 1 Ol [2 1 2 ; 
Mivaii et al]|201 U lAlievato et al.ll201 ll 120121: iKoutoulidis et alj 
2013' lMountrichas et al.ll2013l) . With detailed knowledge of the 
behavior of the linear bias factor with the redshift and luminosity 
we will be able to improve our understanding of major questions, 
such as what is the environment of AGN, what are the major 
triggering mechanisms of AGN activity and how super massive 
black holes (SMBH) co-evolve with the DMH over cosmic time. 



3.1. Method 

The linear bias factor is measured by comparing the amplitudes 
of the observed PS of tracer objects and of the theoretical PS 
of the DM, under the assumptions of a certain cosmology. Since 
the PS amplitude of tracer objects is proportional to the square of 
the linear bias factor (oc b 2 ), its uncertainty directly reflects the 
uncertainty of measuring the latter. Knowing the amplitude (A) 
of our angular power spectrum we are able to estimate this un- 
certainty. The S/N for measuring the normalization of the power 
spectrum Q assuming that its shape is known, is given by: 



- - A 
~N ~ ~5A 



(z) / n \ 2 
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Fig. 3. Signal-to-noise ratio of the amplitude of the angular 
power spectrum (Eq.|5j as a function of the redshift for different 
sky fractions. A Az — 0.2 binning is assumed. 
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Fig. 4. Signal-to-noise ratio of the amplitude of the angular 
power spectrum as a function of the survey duration for the red- 
shift bin 0.8 < z < 1.0 at different sky fractions. 



Here, we are assuming that all multipoles are independent. 

We are using a redshift binning of Az = 0.20 for our cal- 
culation (see Fig. [3] and [5]) but other bin sizes would also be 
applicable to demonstrate our results. In current observations 




(5) 
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Fig. 5. Same as Fig.[3]but for different luminosity ranges (in units 
of log(L[erg s~'])) and the sky coverage of 10 000 deg 2 . 



3 



A. Kolodzig et al. (20!3): AGN and QSOs in the eRASS 



(e.g. [Allevato et al.ll20 1 li iKrumpe et ai]l20 1 2t iKoutoulidis et aD 



2013) a typically much larger bin size is used in order to achieve 
a reasonable S/N ratio for b in each redshift bin. 

3.2. Results 

In Fig. [3] the achievable S/N of the power spectrum amplitude 
is shown as a function of redshift. The shape of the curves is 
dominated by the redshift distribution of AGN modified by the 
quadratic-like increase o f the linear bias fa ctor with redshift at 
constant DMH mass (e.g. lSheth et alj |2001). We can see that we 
are able to measure the amplitude to a high accuracy (< 10 %) 
for a wide redshift range even with a rather small fraction of the 
sky (e.g. ~ 2 500deg 2 ). 

The analysis of the linear bias factor can be performed before 
the entire four year long eRASS is completed, as we can see 
from Fig. g] For a SDSS-like sky coverage of 10 000 deg 2 (blue 
curve) one can already work with the data of eRASS after only 
1 .5 years (three full sky scans) to study the evolution of the linear 
bias factor to an accuracy of better than 10 % in the amplitude for 
the redshift bin 0.8 < z < 1.0. For a sky region of 2 500 deg 2 it 
needs five full sky scans (2.5 years). For the neighboring redshift 
bins 0.6 < z < 0.8 and 1.0 < z < 1.2 the results are similar. The 
sensitivities used fo r this calculation are taken from Fig. 1 of 
lKolodzigetaH(l2012h . 

Owing to the high S/N of the power spectrum amplitude, 
we will be able to separate the AGN in different luminosity 
groups. This is demonstrated in Fig. [5] for a sky coverage of 
10 000 deg 2 . We can see that we will be able to achieve an accu- 
racy of < 10 % for most luminosity groups for a certain redshift 
range. This means that it will be possible to perform redshift 
and luminosity resolved analysis of the linear bias factor with 
eRASS with high statistical accuracy. We note that in our cal- 
culation the difference in the S/N of the luminosity groups in 
Fig.|5]is driven only by the difference in the redshift distribution. 
Although there are some observational evidence for a correlation 
between the DMH mass and the luminosity of AGN/quasars (e.g . 
Krumpe et aT] 120121: iRichardson et all l2012t IShen et alJ 120121: 
Koutoulidis et alJ[2013h . it is still rather tentative. We therefore 
assumed that the source luminosity is independent on the DMH 
mass. 



4. Baryonic acoustic oscillations 

Acoustic peaks in the power spectra of matter and CMB radia- 
tion are one of the major probes to measure the kinematics of the 
Universe (e.g. IWemberg et ai1l2012|). Thev were predicted the- 
oretically over four decades ago (ISunvaev & Zeldovichlll970t 
lPeebles& Yulll970h and now have become a standard tool of 
observational cosmology. Unlike acoustic peaks in the angular 
power spectrum of CMB, their amplitude in the matter power 
spectrum in the ACDM Universe is small. For this reason, 
only recently galaxy surveys have reached sufficient breadth 
and depth for the first convincing dete c tion of BAO, achieved 
with the SPSS data JCole et al. 120051: lEisenstein etaD l2005t 
lHutsill2006t iTegmark et al .1120061) . Since then BAO have been 
measured extensively up to the redshift z ~ .8, in particular 
with luminous red giant galaxies (LRGs) (e.g. An derson et alJ 
[20T2h . Above this redshift limit, BAO features were only found 
in the correlation function of the transmitt ed flux fraction in 
the Lyman-q fore st of high-redshift quasars (Busca et al. 2013; 
ISlosar et al.l2013l) . but have not yet been directly detected in dis- 
tribution of galaxies. 
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Fig. 6. Baryonic acoustic oscillations in the power spectrum for 
the extragalactic sky in redshift range of 0.0 < z < 3.0. At 
wavenumbers above the vertical dashed line (corresponding to 
0.2 h Mpc~') our assumption of a linear clustering starts to break 
down. The red curve shows the original input model for the 
BAOs. 



For the so far uncharted redshift range from z ~ 0.8 up to 
~ 2.0, AGN, quasars and emission-line galaxies (ELGs) are pro- 
posed to be the best tracers to measure BAOs, however, currently 
existing surveys do not to achieve the required statistics for a 
proper detection (Sawangwit et al. 2012; Comparat e t al.ll2013l 
eRASS and the proposed SDSS-IV0 survey program eBOS 
(2014-2020) will be the first surveys to change this situation. 
eRASS will achieve sufficiently high density of objects N ~ 
40 deg~ 2 in this redshift range and will have by far the largest 
sky coverage in comparison to eBOSS and all other proposed 
BAO surveys. This would enable one to push the redshift limit 
of BAO detections in the angular power spectra of galaxies far 
beyond the present limit of z ~ 0.8. 

4.1. Method 

By construction of t he model, BAOs a re included in the AGN 
clustering model of iHutsi et alJ d2012l) . through the 3D linear 
power spectrum (Sect. [2). Oscillations can be noticed, for exam- 
ple, in Fig. [2] in the power spectra of objects selected in narrow 
redshift intervals. As the angular scales of acoustic peaks depend 
on the redshift, in the power spectra computed for broad redshift 
intervals BAO are smoothed out, due to superposition of signals 
coming from many different redshift slices. 

Although analysis of the real data will be conducted in a 
much more elaborate way, for the purpose of this calculation we 
will use a simple method to estimate the amplitude and statisti- 
cal significance of the BAO signal detection. We divide a broad 
redshift interval into narrow slices of the width Az and for each 
slice compute the angular power spectrum, Q(z), and convert 
multipole number to wavenumber k = £/r{z) to obtain P(k,z). 
These power spectra are co-added in the wavenumber space to 
obtain the total power spectrum P(k) of objects in the broad red- 
shift interval. The so constructed power spectrum will have un- 
smeared BAO features. To estimate their statistical significance, 
we also construct a model Cf, smoot h(z) without acoustic peaks, 
by smoothi ng the matter transfer f unction, similar to how it 
was done in Eise nstein & Hul <[[998). From this model we com- 



5 http://www.sdss3.org/future/ 

6 http : //lamwws . oamp . f r/cosmowiki/Pro j ect_eBoss 
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pute the smoothed power spectrum in the wavenumber space, 
^smoothW, not containing BAOs. To illustrate the amplitude of 
the BAO signal, one can plot the difference P(k) - P smoo th(k) or 
the ratio P(k)/P smoot h(k). 

By the analogy with Eq. ©, the S/N of the BAO detection 
in the eRASS sample can be computed as: 



S_ 
N 



/ Q(z) - Q,smooth(z) ^ 2 



(6) 



where the outer summation is performed over the redshift slices 
and the variance cr^, is calculated from Eq. [4] 

The result of this calculation depends on the choice of the 
thickness of the redshift slice Az. For too large values of the Az 
BAO will be smeared out, as discussed above (cf. Fig. |2j. On the 
other hand, for too small values of the Az, at which the thick- 
ness of the redshift slice becomes somewhat smaller than the 
comoving linear scale of acoustic oscillations, the cross-spectra 
between different redshift slices will need to be taken into ac- 
count in computing the P(k). For the purpose of these calcula- 
tions we chose Az = 0.05. The corresponding thickness of the 
redshift slice at z ~ 1 is approximately equal to the comoving 
linear scale of the first BAO peak. Note that the omission of the 
cross-spectra in our calculation leads to some underestimation 
of the confidence level of BAO detection. 

4.2. Results 

In Fig. [6] and [7] we show the ratio P(k)/P smoot h(k) along with 
its uncertainties computed as described above. As these plots 
demonstrate, with the whole eRASS AGN sample for the extra- 
galactic sky we should be able to detect the BAOs with a con- 
fidence level (CL) of ~ 10cr (Fig. |6}. For the currently unex- 
plored redshift range of 0.8 - 2.0 the confidence level of ~ 8<x 
will be achieved, which can be seen in the top panel of Fig. [7] 
Decreasing the sky area to 20 000 deg 2 or 10 000 deg 2 (see mid- 
dle panel of Fig. |7]l we obtain ~ 6cr and ~ 4<x, respectively. In 
Fig. |8]we show for different redshift ranges how the confidence 
level of the BAO detection depends on the sky coverage. The 
curves follow a / sk y' 5 - dependence, as it is expected from Eq. @. 

As Fig.[8]shows, for the redshift ranges 0.0 - 0.8, 0.8 - 1.2 
and 1.2 — 2.0 the achievable confidence levels are rather similar, 
therefore the power spectra ratio shown in the bottom panel of 
Fig- His representative for all three redshift intervals. Comparing 
the upper and bottom panels of Fig. [7] one can see how the BAO 
signal depends on the redshift range, whereas the top and middle 
panels show the degradation due to reduced sky coverage. 

5. Discussion and conclusions 

5.1. Linear bias factor 

Measurement of the linear bias factor provides a simple and di- 
rect method to estimate the average mass of DMHs hosting a 
given sub-population of AGN. With the eRASS, these measure- 
ments will become possible to unprecedented detail. The dra- 
matic improvement of the redshift and luminosity resolution of 
DMH mass measurements will have a great impact on our un- 
derstanding of the environment of AGN, AGN triggering mech- 
anisms and SMBH co-evolution with the DMH. 

An observational evidence of a correlation between DMH 
mass and the luminosity of AGN has already been found but 
uncertainties are still large and AGN luminosities available for 
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Fig. 7. Same as Fig. [6] but for different redshift ranges and sky 
coverages. 
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Fig. 8. The confidence level of a BAO detection as a function 
of sky coverage for different redshift ranges (see Sect. 14.11 for 
further explanations). The vertical gray dashed line shows the 
area of the extragalactic sky. 
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pared their results from clustering studies of AGN in four extra- 
galactic X-ray surveys of different depth and coverage (CDFN, 
CDFS, COSMOS and AEGIS) with theoretical predictions of 
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iFanidakis et al.l d2012l) . Their goal was to determine, from the 
AGN bias factor measurements, the dominant SMBH growth 
mode for AGN of different luminosities - either through galaxy 
mergers and/or disk instabilities or through accretion of hot gas 
from the halo of the galaxy. However, uncertainties of bias fac- 
tor measurements and of DMH mass estimation were too large 
to clearly disentangle the dominant growth mode as a func- 
tion of luminosity. In particular, the luminosity range of ob- 
jects available for their analysis, L ~ 10 42 ~ 44 erg s" 1 , was too 
narrow to challenge the prediction of Fanidakis et al] (1201 2l) 
that luminous galaxies with L > 10 44 erg s _I reside in DMHs 
of moderate mass of ~ 10 I2 M o . For the same reason, no 
direct comparison was possible with the resul ts of the opti- 
cal quasar surveys (Alexander & Hickoxll20l2l) . To overcome 
this limitation, lAllevato et al.l ( 201 1 ) studied broad line (BL) 
AGN from the COSMOS survey and found for their luminos- 
ity bin L » 10 43 ~ 46 erg s _1 a significantly higher DMH mass 
than inferred from quasar studies, suggesting that for broad 
line AGN major merger may not be the dominant triggering 
mechanism, in reasonable agreement with recent s imulations 
dDraper & Ball antvne 2 0T2t iHirschmanneT al. 2012). However, 
the large width of the luminosity bin required to accumulate suf- 
ficient statistics did not allow them to draw a firm conclusion. 

As illustrated by Fig. [5] the eRASS AGN sample will not 
only dramatically improve statistics but will also expand the lu- 
minosity range beyond L ~ 10 44 erg s _1 , to the luminosity do- 
main characteristic of quasars. Thus, the eRASS data will not 
only increase the redshift and luminosity resolution of DMH 
mass estimations of AGN, but will open possibilities for detailed 
comparison of clustering properties of luminous AGN and opti- 
cal quasars. Another aspect of bias measurements with eRASS, 
determining their uniqueness is that they are based on the X-ray 
selected AGN sample and will cover a very broad SMBH mass 
range, broader than t hat in AGN sample s produced by optical/IR 
or radio surveys ("e.g. lHickox et al.ll2009l) . 

The growth rate of SM BHs over ti me can be measured from 
the XLF of AGN (e.g. lAirdet al.l2010l) and eRASS will improve 
the accuracy and redshift resolut ion of these studies tremen- 
dously (e.g. iKolodzig et alj l2012f) . Combined with clustering 
bias data, these measurements will be placed in a broader con- 
text and connected with DMH properties, and will thus provide 
new insights on the co-evolution of SMBHs with their DMHs 
and will also help to investigate the dominant triggering mecha- 
nisms of AGN activity. 

The AGN clustering model used in this paper and, corre- 
spondingly, calculations of the AGN linear bias factor ignored 
the internal structure of DMHs, i.e. they were restricted to the 
scales larger than the size of a typical DMH. Expressed in the 
language of the halo occupation distribution (HOD) formalism, 
these calculations operated with the population averaged halo 
occupation numbers. The angular resolution of the eRosita tele- 
scope, ~ 30" FOV averaged HEW, is sufficient to resolve sub- 
halo linear scales. Clustering measurements on small scales will 
permit to obtain a detailed picture of how AGN are distributed 
within a DMH (e.g. to measure fractions of central and satellite 
AGN) and how the HOD depend on the DMH mass and redshift, 
and AGN lumino isty. Extrapolating result s of XMM-COSMOS 
data analysis by iRichardson et al.1 d2013l) we may expect that 
high accuracy determination of the HOD parameters will be eas- 
ily achieved with eRASS data, which will be able to address all 
these questions, advancing our understanding of AGN clustering 
on small scales and their HOD. 
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Fig. 9. The effective volumes of BAO surveys listed in TableQ]as 
a function of the wavenumber. Effective volumes are computed 
for redshift ranges indicated in the plot. 



5.2. BAO 

The BAO detection beyond redshift ~ 0.8 will be a very signif- 
icant milestone for the direct measurement of the kinematics of 
the Universe. eRASS will be able map this uncharted redshift 
region with a sufficiently high AGN number density to mea- 
sure BAOs with a high statistical significance (see Fig. [8). For a 
proper forecast of how these measurements would improve our 
constraints on cosmological para meters, Markov Chain Monte 
Carlo (MCMC) simulations (e. g. iLewis & Bridld l2002h and/or 
Fisher matrix calculations (e.g. iTegmark et al.lll997|) have to be 
made, which is beyond the focus of our work. Sawangwit et al. 
(1201 2l) have performed a MCMC simulation for QSOs and 
demonstrate that a 3 - Act BAO detection (of a 3 000 deg 2 QSO 
survey with N — 80 deg -2 ) for 1 .0 < z < 2.2 can already signifi- 
cantly reduce the uncertainties. Although the survey parameters 
of eRASS differ (much larger sky coverage but smaller source 
de nsity for the same red shift region, ~ 40 deg 2 ), the results 
of Sawan gwit et al.l d2012l) can give one an idea of how eRASS 
AGN sample will improve the accuracy of cosmological param- 
eters determination. 

Our calculation are limited to the linear regime, and do 
not take non-linear structure growth into account, which would 
somewhat smear out the BAO signal. T his will lead to a decreas e 
of their detection significance (e.g. lEisenstein et al.l 120071) . 
However, with BAO reconstruction metho ds one will be able 
to correct for this effect to s ome extent (e.g. Padmanabhan et aD 
20 12; I Anderson et al.l2012l) . We should also note, that our confi- 
dence level estimation of the BAO detection are rather conserva- 
tive as they neglects information contained in the cross-spectra. 
This will counter balance the negative effect of BAO smearing, 
as it will be d emonstrated in a forthcoming paper dHiitsi et al.l 
12013. in prep.l) where we will study BAO forecasts in a broader 
context and with accurate account of these effects. 



5.2.1. Comparison with BAO surveys 

We now compare the potential of the eRASS AGN sample 
with dedicated BAO surveys in the optical band. For the lat- 
ter, we consider the a lready completed BOSS CMASS survey 
(lAnderson et al.ll2012l) . the planned eBOSS survey which im- 
plementation starts in 2 014 and will continu e till 2020 and the 
future BigBOSS survey dSchlegel et al.ll201 ll) anticipated in the 



6 



A. Kolodzig et al. (20!3): AGN and QSOs in the eRASS 

Table 1. Parameters of BAO surveys. 



Survey 


Tracer 


btr 


Redshift 


^survey 


N 


(n) 


V eff [h" 3 Gpc 3 ] 


Implem. 


Ref. 




Object 


(z = 0) 


Range 


[10 3 deg 2 ] 


[deg- 2 ] 


[10- 4 h 3 Mpc- 3 ] 


(k = 0.07/i Mpc" 1 ) 


Date 




eRASS 


AGN 


- 1.33 


0.8-2.0 


~ 34.1 


-40 


- 0.12 


-7.8 


2014-2018 




BOSS 


LRG 


~ 2.00 


0.4-0.7 


~ 10.0 


~ 80 


- 2.3 


- 2.9 


finished 


[1] 


eBOSS 


ELG 


~ 1.00 


0.6-1.0 


~ 1.5 


~ 180 


- 2.1 


- 0.5 


2014-2020 


|2] 


eBOSS 


QSO 


~ 1.20 


1.0-2.2 


~ 7.5 


-90 


- 0.21 


- 1.7 


2014-2020 


|2] 


BigBOSS 


ELG 


~ 0.84 


0.7-1.7 


~ 14.0 


- 1730 


-6.3 


- 24.0 


> 2020 


[3] 



Notes. b tI - bias factor of the tracer object; f2 slllvcy - solid angle covered by the survey N - surface number density; {n) - average volume number 

density; V cS - effective survey v olu me at the first BAO pe ak 

Ref: ril lAnderson et al] <2012h and lDawson et al.1 d2013h : [2] eBOSS team, priv. comm.; m lSchleeel et al]d201lh 



2020 timeframe. The Table Q] summarizes key parameters of 
these surveys relevant for the BAO studies. 

A quantity often used to estimate the statistical perfor- 
mance of a galaxy c lustering survey is its effective volume (e.g. 
Eisenstei net alJl2005h : 



VeffW = a, 
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n(z)P a (k,z) 



n(z)P tr (k,z) + l 



dV(z) 
dzdQ 



dz (7) 



where n surV ey is the solid angle covered by the survey, P tr (k, z) = 
btr(z) (g(z)/g(Q)) P(k) is the power spectrum of objects used as 
LSS tracer, b ti (z) is their redshift-dependent bias factor and n(z) 
is their redshift distribution [h 3 Mpc -3 ]. Other quantities are de- 
fined in the context of Eqs. (Q~|)-([3]). For optical surveys we will 
assume that b a (z)g(z) = constant, therefore we only need to 
compute P tI (k, 0). The n{z) dependences for optical surveys were 
taken from references listed in Table [1] 

The results of these calculations are plotted in Fig. [9] where 
we show the effective volumes of different surveys as a func- 
tion of the wavenumber. Their values at the first BAO peak are 
listed in the respective column of Table Q] In these calculations 
integration in Eq. (0 was performed over the optimal redshift 
range of each survey, as listed in Table [1] For eRASS we used 
z = 0.8 - 2.0 range in order to emphasize its strength in this un- 
charted redshift region. The eRASS effective volume for the full 
redshift range is ~ 70% larger for the first BAO peak. 

The result of effective volume calculations obviously de- 
pends on the assumptions regarding values and redshift depen- 
dences of comoving density, bias and the growth factor, which 
are not always precisely known, especially for the future sur- 
veys. Nevertheless, these curves should give a reasonably accu- 
rate comparison of how good different surveys are in measuring 
the power spectrum at different scales (note that the uncertainty 
of the power spectrum is proportional to V^' 5 )- 

We can see from Fig.|9]that the effective volumes of eRASS 
AGN and eBOSS QSO samples fall most rapidly towards 
smaller scales in comparison with other surveys. This is a conse- 
quence of the lower volume density of X-ray selected AGN and 
optical QSOs (Table[TJ. For the same reason, the statistical errors 
in the eRASS AGN power spectrum are dominated by the shot 
noise, but the high sky coverage of eRASS keeps them small. As 
one can see from the figure, eRASS is more competitive at larger 
scales, upto the second BAO peak, where its sensitivity becomes 
comparable to BOSS. It should be noted however, that the BOSS 
AGN sample covers a relatively low redshift domain, z 5= 0.7, 
whereas all other surveys presented in Fig. |9]are aimed at sig- 
nificantly larger redshifts, z ^ 0.7 (Table [TJ. Around the first 
peak the effective volume of eRASS is a factor of ~ 2 - 3 larger 



than for BOSS but a few times smaller than BigBOSS (TablefB- 
On the other hand, eRASS exceeds eBOSS at all wavenumbers. 
This would still be the case if one considers subset of the eRASS 
sample covering only - 1/3 of the extragalactic sky. 

In the conclusion of this section we note that it is remarkable 
that the statistical strength of eRASS for BAO studies is compet- 
itive with that of dedicated BAO surveys, even though eRASS 
was never designed for this purpose. Potentially, the eRASS 
AGN sample will become the best sample for BAO studies be- 
yond redshift z <: 0.8 until arrival of BigBOSS in the end of 
this decade. However, this potential will not be realized without 
comprehensive redshift measurements. 

5.3. Redshift data 

We assumed so far that redshifts of all eRASS AGN are known. 
Now we briefly outline requirements to the redshift data imposed 
by the science topics discussed above. 

The linear bias as well as luminosity function studies do 
not demand high accuracy of the redshift determination. Indeed, 
values of the order of & - 0.1- 0.2 should be sufficient, un- 
less analysis with much higher redshift resolution is required. 
In principle, such accuracy can be provided by photometric sur- 
veys. However one would need to investigate the impact of large 
fraction of catastrophic errors, from which AGN redshift deter- 
mination based on the standa rd photometric filter sets are known 
to suffer (Salv ato et al.ll2.01 Th . Of particular importance are red- 
shift and luminosity trends in catastrophic errors. These prob- 
lems will be c onsidered in the forthcoming paper (Htits i et al] 
12013. in prep.l) . Provided that they are properly addressed, opti- 
cal photometric surve ys of a moderate depth of / ^ 22.5 mag 
dKolodzig et al] 120121) and with the sky coverage exceeding 
2 500 deg 2 (Fig. [3]l would already produce first significant re- 
sults. An existing survey with such parameters is SDSS. Its depth 
would allow detection of w 80% of eRASS AGN (Kolodzi g et al] 
2012) and with its sky coverage of - 14 500 deg 2 one should 
be able to conduct high accuracy measurements of the linear 
bias factor. A mong other, on-going surveys, the Pan-STA RRS 
PS1 3tt survey dChambers & the Pan-STARRS Teamll2006l) ful- 
fills the necessary depth and sky coverage criteria. 

BAO studies, on the other hand, require a much higher 
redshift accuracy of the order of dz ~ 0.01. Such an accu- 
racy can be only achieved in spectroscopic surveys or in high 
quality narrow-band multi-filter photometric surveys. For exam- 
ple, for a 4<r detection of BAOs in the redshift range 0.8 < 
z < 1.2, a spe ctroscopical survey is needed of the depth of 
/ > 22.5 mag (Kolodzi g et al] 120121) and sky coverage of at 
least ~ 20 000 deg 2 (Pig. [8J. Promising candidates are the pro- 
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posed 4MOST (Ide Jong et alJl2012t) and WEAVE dDalton et alJ 
120 121) surveys, which would cover a large fraction of the sky 
with a multi-object spectrograph in the southern and northern 
hemisphere respectively. An important caveat is that the angular 
resolution of eRASS (FOV averaged HEW of ~ 30") is insuf- 
ficient to provide accurate positions of X-ray sources for spec- 
troscopic follow-ups with multi-object spectrographs. Therefore 
additional photometric surveys (for example Pan-STARRS PS 1 
3n) will be needed in order refine source locations to the required 
accuracy. 

Finally, we note that we excluded from consideration a num- 
ber of observational effects and factors, such as source confusion 
and source detection incompleteness, positional accuracy, tele- 
scope vignetting and non-uniformities in the survey exposure, 
and a number of others. These factors and effects are well known 
in X-ray astronomy and data analysis methods and techniques 
exist to properly address them in the course of data reduction. 



6. Summary 

We have explored the potential of eRosita all-sky survey for 
large scale structure studies and have shown that eRASS with 
its « 3 million AGN sample will supply us with outstanding op- 
portunities for deta iled LSS research. O ur results are based on 
our previous work dKolodzig et al.ll2012l) . where we investigated 
statistical properties of AGN in eRASS, and the AGN clustering 
model of iHiitsi et al] d2012h . 

We have demonstrated that with eRASS, the linear bias fac- 
tor of AGN can be studied to unprecedented accuracy and detail. 
Its redshift evolution can be investigated with an accuracy of bet- 
ter than ~ 10 % using data from the sky patches of ~ 2 500 deg 2 . 
Using the data from sky area of 10 000 deg 2 , statistically ac- 
curate redshift and luminosity resolved studies will become pos- 
sible for the first time. Bias factor studies will yield meaningful 
results much before the full 4 years survey will be completed. 
The eRASS AGN sample will not only improve the redshift and 
luminosity resolution of bias studies but will also expand their 
luminosity range beyond L0.5-2.0keV ~ 10 44 erg s _1 , thus mak- 
ing possible direct comparison of clustering properties of lumi- 
nous X-ray AGN and optical quasars. These studies will dra- 
matically improve our understanding of the AGN environment, 
triggering mechanisms, growth of super-massive black holes and 
their co-evolution with dark matter halos The photometric red- 
shift accuracy should be sufficient for the bias factor studies, al- 
though the impact of large fraction of catastrophic errors typical 
for standard broad band fi lter sets yet needs to be investigated 
dHiitsi et alJl2013. in prep]) . 

For the first time with the X-ray selected AGN, eRASS will 
be able to detect BAO with high statistical significance of ~ 10cr. 
Moreover, it will push the redshift limit of BAO detections far 
beyond the current limit of z ~ 0.8. The accuracy of BAO in- 
vestigation in this uncharted redshift range will exceed that to 
be achieved by eBOSS, planned in the same timeframe and will 
be only superseded by BigBOSS proposed for implementation 
beyond 2020. Until then, eRASS AGN can potentially become 
the best sample for BAO studies beyond z ^ 0.8. However, in or- 
der for this potential to be realized and exploited, spectroscopic 
quality redshifts for large fraction of the sky will be required. 
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